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ABSTRACT 

The transient testing technique for the determination of heat trans- 
fer characteristics of a porous medium has been known since 1929, It 
requires a good deal of care and one of the important restrictions is a 
step change in temperature must be applied to the medium at the beginning 
of a test run. This thesis involves the design, test and evaluation of 
an experimental set-up that utilizes an "integral heater" to accomplish 
the step change in temperature. In addition, there is presented heat 
transfer and flow friction characteristics for a 20° skew matrix con- 
structed from perforated nickel material. This data is compared to a 
20° skew matrix made from stainless steel material that has been pre- 
viously reported by C. P. Howard. Also, data is presented for two "paral- 
lel plate" matrices, one constructed from brass sheet, the other from 
perforated nickel material. A comparison of their relative performances 


is made, 
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Flow cross-section area, ft 
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Matrix mass, 1b 
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Fluid viscosity, lb/hr ft 
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Matrix density, 1b/£t> 
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Fanning friction factor, ratio of the wall shear stress 
to the fluid dynamic head 
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Colburn j factor, Noor 
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Reynold's number 4x. G/ 


Stanton number, see 
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Number of heat transfer units, hA/ Tic, 
Entrance flow coefficient for the matrix 
Exit flow coefficient for the matrix 
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Matrix flow void volume/Matrix volume 


Conduction parameter, kA, / Linc, 


Subscripts: 
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upstream 
downstream 


fluid 


vili 





solid (matrix material) 


Standard temperature and pressure 





le INDRODUCTION, 

A compact heat exchanger is one not necessarily of small bulk and 
weight, but one which incorporates a heat transfer surface having a high 
area density. 

Arbitrarily, it has been specified that & = 200 Fone makes 
for a compact heat exchanger, where B is the ratio of the heat trans- 


aN 


fer area to core volume, ft’ /ft 
Motivation for using compact surfaces is to gain high exchanger 
effectiveness in given box volume and weight limitations and at a reason- 
able cost. Great impetus has been given the compact heat exchanger field 
by recent applications as a regenerator in aircraft turboprop power plants, 
as a regenerator in the gas turbine automobile engine and by their use 
in cryogenic plant exchangers in the production of liquid oxygen and 
nitrogen, As the applications for compact heat exchangers increase, 
design problems increase, A successful compact heat exchanger design 
and development program requires access to, or the ability to generate 
knowledge of the basic heat transfer and friction behavior of compact 
surfaces of interest, It is also desirable that the designer would have 
available as much information on surface performances as possible to re- 
main flexible in his design with regard to the limitations that have been 
placed upon him, 
Building a complete heat exchanger to evaluate a new surface as to 
its heat transfer and friction performance would be prohibitively ex- 
pensive. The alternate method is to generate the performance character- 


istics of a surface of interest by testing small matrix modules and use 


+Numbers in parentheses refer to references in the bibliography. 
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this data to design a prototype heat exchanger. 

Briefly, the objectives of this thesis are: to build an experimental 
rig that will test small matrix modules to determine their heat transfer 
and flow friction characteristics, and, to test matrices made from per- 
forated nickel material to determine these characteristics. 

One of the simplest and most satisfactory methods of measuring the 
forced convection heat transfer characteristics of a matrix utilizes a 
transient test technique or so called single blow method. This test 
technique consists of heating the homogeneous porous solid of known 
weight and geometry with a uniformly distributed heat source and then 
at a given instant, removing the heat source and simultaneously cooling 
the solid with a fluid at a uniform lower temperature. The heating and 
cooling steps can be reversed and the same results achieved. A time- 
temperature response curve of the fluid stream leaving the solid is 
recorded, Ry comparing this record to theoretically developed non-dimen- 
sional temperature response curves, the heat transfer characteristics of 


the solid can be determined. 





Be SUMMARY OF THEORY. 

A theoretical analysis of the transient heat transfer behavior of 
a porous medium was presented by Schumann (2) in 1929. <A porous medium 
may be defined as an aggregation of solid bodies so arranged that the 
void spaces between the bodies are connected and will permit the passage 
of a fluid. This definition can of course be extended to include the 
matrix surfaces utilized in a compact heat exchanger. 

In Schumann's analysis, the porous medium is considered to be homo- 
geneous and initially at a uniform temperature, ti. A fluid at the same 
temperature is flowing through the medium, At a certain instant the 
temperature of the entering fluid is assumed to change to a higher tem- 


perature t The problem is then to find the temperature of the fluid 


oe 
and the solid as functions of time and position in the matrix or porous 


media. The basis of the analysis are the theoretical generalized heat- 


ing curves that are developed from an energy balance on the matrix. 





Bek, ee hb(t,-te)dx age rep (t,- Sdx) 
: at, 2’t. 
be — Pua be ba Ce A 





1) Properties of the fluid are independent of temperature. 

2) Fluid flow is steady. 

3) The porous solid is homogeneous. 

4) The thermal conductivity of both fluid and solid is infinite 
perpendicular to the flow direction. 

5) Thermal conductivity of fluid is zero in the flow direction, 


ats 
Heat absorbed by the solid = Ps As CG. 96) ax 
Heat transferred to solid by convection = na 
Heat transferred from the fluid by convection = MCo a oor 


: 
Heat transferred from solid by conduction = =~ k, Ne <4 dx 


An energy balance for the fluid and solid results in: 
aes 
p,As C, ae Se dx = lie On ok + hb (4, - ts) dx (1) 


Mm Cp ae dx + hb(t,-ts)dx (2) 


Oo 
Let: T= generalized time variable bA Ss) 
Ss 





Ss 
generalized position variable = a 
tu L 
N = h A De 
tu m c Cisse 
P m E L 


Making the above substitutions, equations 1) and 2) become; 
ats ae 
ab = (4, —ts)+ ANeu Sas 
= (t,- t+) 


The experimental method used requires that the above assumptions 





and the following boundary and initial conditions be closely observed. 
These conditions are as follows: 
1) Initially the core is at a uniform temperature. 


2) At time t = 0, the temperature of the incoming fluid changes 
instantaneously. 





sll 


3) System is adiabatic, that is, only heat is transferred be- 
tween the matrix core and the fluid,none with the surroundings. 


With the idealizations, initial and boundary conditions, and the additional 
assumption that the thermal conduction in the solid is zero parallel to 
the direction of fluid flow (i.e., A , the material conduction para- 


meter, = 0), Schumann's solutions are as follows: 


[> @) 

(en = fe n at 

ee ees 2 acary") «Je (ni V2)... 
an 


fl “i 
CO 
t. Sate 3 
eee (2 +7) 22 —2___ 3, (2i 2) 
fie ot d(z T ) 
n=O 


This solution was first used as the basis for a transient technique by 
Furnas (16) in 1932, The fluid temperature was measured at the exit of 
the matrix where x = L, therefore z = Nee The experimental data was 
compared with Schumann's theoretical constant z curves. The theoretical 
z curve which best fit the experimental data was considered as the ee 
of the matrix and from this value, h, the convection conductance was 
determined, The biggest disadvantage of this method is in the calcula- 
tion of an infinite series for each point and the number of terms neces- 
sary for a good approximation increases as Sie increases, 

In 1950, G. L. Locke (1) developed another method for the evaluation 
of experimental He from the theoretical curves. In this method, only 
the maximum slopes of the theoretical and experimental curves need to be 
compared, He derived an expression for the slope of the generalized heat- 


- - 1 U 
ing curve, (t, t / (te, t.) vs. T/z, as a function of Neu Locke's 





expression shows that the maximum slope of the generalized heating curve 
is a unique function of Ne 


a 


4 2 
deren aaeds) Neu 
t 


t = (. = ; -(N cep 
£1 i = IN rt - J,(2i Irn ) e tu 
d(T IN, tu 1 tu 


os 


The above implies that the downstream fluid temperature is measured at 
x = L. Therefore, z = Ne and te = tro: This method has become known 
as "The Maximum Slope Method". 

One of the idealizations made by Schumann and followed by Locke was 
that the thermal conductivity of the solid parallel to the flow is zero. 
This is very difficult to achieve in an experiment with a continuous flow 
channel in the matrix. Although the effect of finite thermal conductiv- 
ity parallel to the flow is negligible at high Reynold's numbers, it can- 
not be ignored at low Reynold's numbers. Recent works by Crestwick (3), 
Mondt (4), and Howard (5) consider the effects of longitudinal thermal 
conduction in the solution of Schumann's transient heating problem. To 
correlate results, they utilize the conduction parameter, A. This 
parameter is the ratio of conduction heat transfer to the change in in- 
ternal thermal energy of the flowing fluid for the same temperature differ- 


ence. When A= 0, the maximum slope is a unique function of N as put 


tu 
forth by Locke, however, when AF 0, the maximum slope is a function of 
N. and A. 
tu 
Howard, in his paper (5), included graphs and a table of maximum 
slopes as a function of Ne and A . He derived these values by finite 


difference calculations on a digital computer and proposed that the values 


of maximum slopes, as presented, are within + 2% for a given phe and A. 
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Ideally then, if the idealizations, initial and boundary conditions 
can be adhered to, all that is required for the experimental evaluation 
of the heat transfer characteristics of a matrix is to: 

1) Determine the conduction parameter Ae 


2) Experimentally determine the maximum slope of the heating or 
cooling curve. 





she OBJECTIVES, 

In the past, there has been a variety of experiments conceived to 
utilize the principles previously discussed, by approaching the idealiza- 
tions, initial and boundary conditions as closely as possible, Perhaps 
the most difficult condition to achieve is; "at time t = 0 the tempera- 
ture of the incoming fluid changes instantaneously". 

Locke (1) and Coppage (9) both employed an air switching device with 
which air could be quickly directed through the matrix or by-passed to 
the atmosphere. By use of this switch they imposed a "step change” in 
temperature on their test matrices. To approach the "step change” in 
temperature, "requires a fast operating fluid switch for introducing 
the hotter (or colder) fluid in order to assure an instantaneous change 
of fluid temperature, However, even with rapid switching, one of the 
major difficulties here is the inevitable thermal capacity of ducting, 
switch device and matrix holding device causes unavoidable changes in the 
entering temperature". (9) 

Solar Aircraft Company recently did some work testing core modules 
for their T-600 engine (10). In their experimental rig, they utilized 
two steady streams of air and a carriage tray that houses the core module. 
The tray was inserted through the side of the test cabinet and was moved 
from a warm air stream to a ambient or cold air stream to simulate a step 
flunction change of temperature. This type of switching was done manu- 
ally. 

Howard, in his work with new compact surfaces (6), employed a "slid- 
ing drawer", This drawer holds the matrix and the downstream thermo- 


couples. The drawer can be inserted or withdrawn from the test section 





through which ambient air flows continuously. The matrix is heated 
when the drawer is in the withdrawn position by an independent tempera- 
ture controlled air stream, 

One of the objectives of this thesis was to construct and test an 
experimental rig that would improve upon those previously used. This 
test apparatus described herein differs from other test facilities 
chiefly in the method used to simulate the initial condition of a step 
change in temperature of the fluid (air). 

The test apparatus constructed has but one air stream, the air is 
heated by a high resistance nichrome wire heater complex, an integral 
part of the apparatus, located just upstream of the test matrix. The 
time constant or response of the heaters is at least an order of magni- 
tude smaller than the response of the thermocouples and recording devices 
used, By turning the heaters on or off, the step change in temperature 
is simulated. The other systems required some manual dexterity on the 
part of the experimenter to begin a data run, whereas this system in- 
volves only turning off or on the heaters. The temperature recording 
device is started prior to the start of a run. The results achieved by 
this system were an improvement over other systems. With the "sliding 
drawer" and "quick changing valve" there is a change in the flow and there- 
fore the hydrodynamic as well as the thermal boundary layer must develore. 
In this system, the hydrodynamic boundary layer is already formed, only 
the thermal boundary layer is affected. A complete description of equip- 
ment is presented in APPENDIX B. 

In addition to the construction of this test apparatus and the verifi- 
cation of its reliability by comparing with published data, extension of 
the curves published by Howard (6) into the high Reynold's number region 


was attempted. 





Finally, heat transfer and flow friction characteristics of several 


new surfaces were determined. 





4, EXPERIMENTAL TECHNICUES 

In reference (5) Howard presents results from a numerical finite- 
difference method of calculating the transient behavior of a porous 
media when subjected to a step change in fluid temperature, These cal- 
culations consider the case where longitudinal thermal heat conduction 
cannot be neglected. These data are intimately involved in the determina- 
tion of ee and subsequently the Colburn j factor as presented in this 
thesis. For this reason the restrictions, idealizations and conditions 
Howard placed on his solution will be restated here. It is these re- 
strictions, idealizations and conditions that the experimental set-up 
must approach as closely as possible to achieve reliable results. 

1) The flow entering the matrix is steady and uniform in velocity 
and temperature, and remains steady and uniform at any cross 
section as it flows through the matrix. 

2) The thermal conductivity of the matrix is infinite in the 
direction normal to the flow and finite in the direction para- 
llel to the flow. 

3) The thermal capacity of the fluid contained at any time within 
the matrix is small compared with the thermal capacity of the 


matrix. 


4) The thermal properties of the fluid and matrix are constant 
and uniform. 


5) The convective heat transfer coefficient is some suitable 
average and remains uniform and constant. 


6) At time zero the change in fluid temperature will be a step 
function, with the matrix and its entrained fluid initially at 
some uniform and equal temperature. 

The first idealization was attained by careful experimentation. 

Uniform velocity profiles were obtained by wire screen flow straighteners. 
Uniform temperature profiles were obtained by a wire screen type heater 


explained in detail in APPENDIX B. The velocity profiles were taken at 


the inlet of the matrix over entire range of flow rates to ensure that 


ll 





the profiles were uniform, The temperature profiles were taken over the 
entire range of flow rates and found to be uniform within + .5°F, 

The second idealization cannot be strictly realized in practice. The 
effect of the deviations from this idealization can be made negligible in 
the solid if the resistance to heat transfer within the solid is small 
compared to the heat transfer resistance between the solid and fluid. 

The third idealization restricts the fluid to a gas rather than a 
liquid. 

The fourth idealization can be realized if the temperature variation 
is limited to 20 - 25°F, rise above ambient, This would place the 
operating temperature range at 60 to 80°F. In this range the thermal 
properties of the fluid, particularly a gas, vary insignificantly and can 
be considered, as a good approximation, constant and uniform. 

The fifth idealization is imposed upon the computer solution and, of 
course, the determination of the convective heat transfer coefficient is 
the ultimate result of the experimental process, 

Finally, two sets of thermocouples were placed one upstream and one 
downstream of the matrix. Prior to the beginning of each run, these 
thermocouple complexes were read. When a steady temperature differential 
was indicated, the run was begun. 

Herein lies the improvement over the "sliding drawer" technique. A 
series of nichrome wire heaters were installed in the air stream to heat 
the air and subsequently the matrix. To begin a run, the heaters were 
turned on or off, and the temperature vs. time recorder started. The time 
constant (TC) of the nichrome wire heater is very small (TC = .0425 sec 


for m = 1000 lb/hr, TC = .259 sec for m = 10 1b/hr) and hence the approach 
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to the final idealization: at time t = 0 the change in Fluid temperature 
will be a step function. 

A schematic diagram of the experimental system in which air is used 
as the working fluid is shown in Figure 1. Figure 2 is a photograph of 
the entire test apparatus. Essentially, the system consists of a source 
of air, an air metering device, flow straightener screen pack, an air 
heater section, a test matrix, a pressure measuring system and a tempera- 
ture measuring system, 

Air is drawn through the apparatus by a six-stage centrifugal air com- 
pressor. Air flow rates are measured by an ASME standard D and D/2 ori- 
fice meter with changeable orifice plates. Mass flow through the system 
is controlled by two sliding valces, one on the inlet to the compressor 
and one on the compressor outlet. In addition, there is a fine control 
by-pass valve located just ahead of the compressor. As the air passes 
through the orifice plate it comes to the screen flow straightener at 
which point the flow cross-section changes from round to square. The air 
then flows through the air heater section which serves a dual purpose of 
straightening the air velocity profile and heating the air. In this 
section the air is heated uniformly to a predetermined temperature, usually 
20°F above ambient. The heated air passes into the matrix and heats the 
matrix to the same uniform temperature. At this point, a run is ready to 
start. By referencing the thermocouples downstream of the matrix to the 
thermocouples upstream of the heaters, the initial temperature differen- 
tial can be very closely controlled. This temperature differential is 
continuously recorded by a Leeds and Northrup ''Speedomax" strip chart 
recorder when a run is started. The recorded trace has a distinct ad- 
vantage in ena there is no transposition of data necessary, hence, a 
decrease in the uncertainty of the results. The resultant trace or curve 
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is the primary data for determining j. The Ny, of the matrix can be evalu- 
ated by determining the maximum slope of this fluid temperature differ- 
ence curve during the cooling transient. This is the basis of the maxi- 
mum slope technique as developed by Locke (1). 

The flow friction performances of the matrices tested were obtained 
by determining the isothermal pressure drop across the matrix and the 
mass flow rate of the fluid through the matrix. The pressure drop data 
were obtained from static pressure taps located in the test section im- , 
mediately upstream and downstream of the test section. The entering 
fluid static pressure was also recorded. The mass flow rate data were 
obtained in the air metering section as previously described. In summary, 
the following data were recorded for each run: 

a) Orifice diameter (D,) 
b) Pressure drop across the orifice ( AP ) 
c) Pressure drop across the matrix ( AP 
d) Entering static pressure (Ps) q 
e) Orifice temperature (T,) 
f) Atmosphere pressure (P.) 
g) Temperature-Time cooling or heating curve 
h) Temperature recorder chart speed 
A complete description of the experimental apparatus and equipment is 


in APPENDIX B. APPENDIX A discusses the relation of the data to the 


equations governing j & f. 
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Se DESCRIPTION OF SURFACES. 

There are data presented in this thesis for seven surfaces. Of 
the seven, three are new. Three of the remaining four surfaces were 
tested by Howard (6). The remaining surfac: ‘was used primarily as a 
test surface to check the reliability of the experimental apparatus. 

The test surface was a matrix consisting of 20 randomly stacked 
stainless steel, 60 x 6C mesh screens. A photograph of the matrix 
along with a tabulation of its physical and geometrical properties may 
be found in Figure 3. The heat transfer and flow friction characteris- 
tics for such a matrix were published by Coppage in 1952 9), and was 
considered a good check on the performance of the experimental set-up. 

The three matrices examined by Howard and re-examined in this thesis 
were: 

1) Plate-fin reference matrix 
2) 10° skew passage matrix 
3) 20° skew passage matrix 
These particular surfaces were chosen to provide an additional check 
on the reliability of the experimental setup and to provide a compari- 
son between the "integral heater" and "sliding drawer" experimental techni- 
ques. 

The plate-fin surface is formed by corrugated sheets separated by a 
flat spacing sheet to prevent meshing of the corrugations, This matrix 
was to serve as a reference, since the theoretical laminar flow solutions 
for both heat transfer and frictional behavior of triangular cross-sec- 
tions are available. However, difficulty was encountered in forming the 
material into triangular cross-sections, and therefore it fails as a re- 


ference, but the results were interesting and will be discussed later, 
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Figure 4 presents photographs of the plate-fin matrix and a tabulation 
of its' physical and geometrical properties. 

The skewed passage also employs corrugated sheets, but unlike the 
plate-fin, the spacer plate is not included. To prevent meshing of the 
corrugations, the sheets are skewed relative to one another forming what 
is defined as the "skew angle". A photograph of the matrix and a tabula- 
tion of its properties may be found in Figure 5. 

Of the three new surfaces tested, one is a version of the previous- 
ly tested 20° skew matrix using perforated or slotted nickel plates for 
a material vice a continuous material. The plates were corrugated on 
the same device used for the original 20° skew to maintain the same corru- 
gation geometry. The perforations should reduce stagnant areas in the 
matrix, thereby increasing the heat transfer characteristics. 

Finally, two matrices of a parallel plate design were constructed 
and tested. The "parallel plate" matrix is similar to the plate-fin 
in that a spacer plate is placed between the formed sheets to prevent mesh- 
ing. The sheets are formed by rolling. When assembled, the matrix is a 
maze of rectangular passages with an aspect ratio of 77. Aspect ratio is 
defined in reference (8) as the ratio of the two sides of a rectnagular 
flow passage, b/a. The sketch below is, ideally, what was intended in 


designing the "parallel plate" matrix. 
os am 
aae/ ee ate 
One matrix was formed from a brass material and is used as a refer- 
ence to theoretical laminar flow data for parallel configurations (8). 
The second matrix was formed from the perforated nickel plate to examine 
the effects of the perforations on this type of geometry. Pictures and a 


tabulation of properties for the "parallel plate" matrices are in Figure 6. 
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6. PRESENTATION OF RESULTS. 
The heat transfer and flow friction data for each matrix are pre- 
sented in Tables I to VII. These data are also presented in graphic 


2/3 


form, the Colburn heat transfer factor, j ( = NowNor }, and the friction 


factor, f, are plotted as functions of Reynold's Number, N in Figures 7 


R? 
through 16. Plotted also in Figures 13 through 15, for comparison pur- 
poses, are the theoretical laminar flow solutions for parallel plates. 
For the heat transfer data, Prandtl number (Nop) was not really a 
test variable, since the working fluid was air at moderate temperatures. 


N zi is included as an approximation to the effect of Prandtl number 


PR 

over a moderate range of Prandtl numbers. Basic laminar boundary layer 

solutions indicate that in the Prandtl number range of 0.5 to 1.0 (gases), 

the Prandtl number enters the solution as approximately the 2/3 power. 
Since a moderate temperature range is used in the experiment, air 

properties throughout the matrix are practically constant. Therefore, 

all the properties of air used in the calculations of the heat transfer 

and friction characteristics were those of air at the bulk temperature of 

the air stream as measured at the orifice inlet. 


The Reynold's number is based on the hydraulic diameter and mass flow 


velocity. 





The hydraulic diameter is defined as; 


D 4 dr Pa Flow Passage Cross-sectional area ) 


H Flow Passage Wetted Perimeter 
= 4 Void Volume x ib: 
L Heat Transfer Area 
aay Cy 
A 


Le 





The mass flow velocity is based on the matrix free flow cross-sectional 


area: 


m m 
A PAr 





The friction factor includes the entrance, exit and flow accelera- 
tion effects, The entrance and exit effects normally provide only a 
smali contribution to the overall pressure drop in a matrix. Core fric- 
tion is the term that controls the magnitude of the pressure drop. As a 
result, high accuracy in the evaluation of the entrance and exit loss co- 
efficients is not required. These coefficients ke and es are functions 
of flow geometry and Reynold's number and were obtained from Chapter IV, 


Figures 18 to 21, as appropriate, reference (8). 


18 





i UNCERTAINTIES AND LIMITATIONS, 

As put forth in the theory, the experimental method utilized, re- 
quires great care to approach the idealizations and initial and boundary 
conditions of the theoretical analysis. The magnitude of errors intro- 
duced if the experimental system does not duplicate the ideal system is 
difficult to evaluate. Each of the conditions and idealizations were 
carefully considered and the experimental system was designed to approach 
as closely as possible the ideal system, These items are discussed in 
the “Experimental Technique" section. No attempt will be made here to 
evaluate or estimate this type of error. 

Other experimental errors are caused by: 1) uncertainty in the 
value of physical constants, 2) inaccuracy in geometric measurements 
and 3) instrument inaccuracies, The method employed herein to describe 
the uncertainty in the values reported is explained fully in Reference (14). 
However, before stating these uncertainties, an example of the process 
follows: 

A theorem of statistics as stated in Reference (14) is: 

"If R is a linear function of n dependent variables, each of which 
is normally distributed, then the relation between the interval for the 


veriables, w;, and the interval for the result, wp, which gives the same 
odds for each of the variables and for the result is: 


Ya 
aR 2 oR 2 aR 2 
We = ce W, ) + (Gu: Wa) focoe fF ay. Wa ) 


where Wee uncertainty interval based on certain odds. 


Consider the Colburn j factor 


7 Ac 2/3 
J= Neu A Npr 
: 2/3 
by the theorem and letting P = Nop 
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. ig 2 * 2 2 oie - 2 
Wis | epee auc + (Mu ALP wy) + (2 Nw Ae 9” wp) 


if both sides are divided by j then: 


= ne Wa. \ * = £ 2 
a Gee) “Gua (4) *(5 ¥) 


Before an equation such as this can be evaluated the uncertainty 


(1) 


= 


terms must be known. 
The physical constants necessary are listed below. References (12) 
and (13) were sources for the values of these constants. The uncertainties 


in these values, as closely as can be determined, are also listed below. 


ec + 0.5% 
S = 

c + 0.5% 
: x 

7 + 1.0% 

Nor + 2.07 

k + 0.5% 
5 é 


The uncertainty in physical linear measurements could be consider- 
ed negligible if core construction did not have minor inconsistencies, 
The uncertainty in the weight of a core is considered negligible. Conser- 


vative estimates of the uncertainty in linear measurements are: 


A, A.» Als Ags + 1.0% 
LL” + 0252 
Wo ale 

Ss 
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Instrumentation errors were such that only the manometers need be 
considered, With the exception of the orifice temperature, all tempera- 
tures measured were temperature differences and these were measured in 
inches as indicated in Appendix A. Since pressures varied and manometers 
to measure them changed, the maximum uncertainty encountered was used in 


the error analysis and is listed below. 


P A257 

5 a 

P + 1.70% 

* = 

Been negligible 


Another uncertainty which requires special mention is that in Neus 


Since j does not vary linearly with N,, the uncertainty at high and low 


R? 
Np was considered. This uncertainty was determined by utilizing Figure 19 


and the uncertainty in Maximum Slope as 2.0%. 


N = 3.5% De 0 oe = + 6.0% 
ee — Fs 
tu 
Ngee 30g A= 004 rane Gerona 
tu ae 
N 
tu 


Returning to Equation (1) and substituting as appropriate: 


. peas x Vo ’ 9 

For high Ne? = [cos (o1)* 4 (. Ol) + sy | & -063 “ BAe 
w. = z z 4 2 Yo : 9.8 7, 

For low Np: a [core +(.01) +01) + cay" | = .098 WwW, - 7. lo 


Following the above type of analysis, the following additional 


uncertainties were calculated, 


m 4.102 
Np 2,07 
£ + 4.3% 


oq 





8. DISCUSSION OF RESULTS. 

The 60 x 60 mesh screen matrix was used solely as a test matrix, a 
matrix whose heat transfer and friction characteristics were well docu- 
mented, Figure 7 and Table I present the results. There is quite good 
agreement in the friction factor but a slight variation in j factor. This 
difference is attributed to the difference in the geometries of the test 
matrices, Exact reproduction of a screen matrix is difficult to achieve. 
The correlation was such that the experimental test setup was considered 
capable of delivering reliable data. 

For further verification on the reliability of the test setup, Figures 
8 to 10 and Tables II to IV present the test results for the Plate-Fin 
Reference Matrix, the 10° Skew Stainless Steel Matrix and the 20° Skew 
Stainless Steel Matrix used by Howard (6). Although this was not the 
original purpose for utilizing these matrices, the results do indicate 
conclusively that the "integral heater" is as good or better than the 
"sliding drawer" technique, 

The original purpose for utilizing these matrices was to try and ex- 
tend Howard's curves into the higher Reynold's number region (lower ee 
region). The points obtained in this region are plotted on the curves and 
the perinent data are recorded in the Tables. None of them follow the 
straight line of their predecessors. There are at least two reasons for 
this obvious departure from the norm. First, it is a well known fact 
that the ''Maximum Slope" Method has an inherently large error in the low 
Ne region of 1.5SN, $3.5. This fact was originally published by Locke 
(1) for A= 0. Howard (5) took into account longitudinal conduction and 


published an error curve that includes the effect of longitudinal conduc- 


tion on an error in New see Figure 19, In the regions below Neu <2 


Ze 





the error is acceptable, yet, a discrepancy still exists. This points 
out the second reason for the discrepancy. The maximum slope at high 
Flow rates (low Ne occurs at time t = 0+. Instrumentation prevents 
accurate measurement of the true maximum slope at this time due to in- 
herent delay or lag time. 

Extrapolation of the time-temperature curves back to time t = 0 was 
attempted in an effort to remove the effects of instrumentation. The 
method used was either incorrect or unreliable. Satisfactory results 
were unobtainable. 

The merits of the plate-fin, 10° skew and 20° skew matrices as heat 
transfer surfaces was discussed thoroughly by Howard in reference (6). 

Figures 11, 13 and 14 and Tables V, VI and VII present the results 
for a new geometry tested and a new idea in material. Figure 11 and 
Table V present the results for a 20° skew nickel matrix constructed from 
slotted plates, Appendix C explains fully the geometry of the slots and 
the calculations that determine the basic geometric parameters p, es ; 
and Ty The effects of longitudinal conduction was considered negligible 
as the matrix was constructed so that the air flow was normal to the long 
dimension of the slots. Plotting all the points uncorrected for conduc- 
tion showed an erratic departure from the expected results for Np < 200. 
It was then decided that perhaps the basic assumption was not warranted 
completely. Corrections for longitudinal conduction were made where 
significant and these points were plotted. These corrected values of j 
are overestimations, for, although A, is corrected for the effects of the 
slots, the matrix length L is not the true conduction path length but at 


best a lesser approximation. Since Xr is obviously not equal to zero, 
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the true value of j is probably between the uncorrected j ( A= 0) and 
the j corrected for longitudinal conduction using L as an approximation 
for conduction path length. This was considered in the interpretation of 
results. 

Figure 12 compares the 20° skew stainless steel matrix with the per- 
forated nickel matrix. There is a 70 to 100% increase in j, while f in- 
creases 10 to 25%. 

It must be understood that these two matrices are geometrically similar 
only in overall dimensions. Actually, the heat transfer areas were cal- 
culated on different bases, The 20° skew stainless steel heat transfer 
area was based on the total surface area available for heat transfer while 
the heat transfer area of the 20° skew perforated was based on this same 
area less the area of the slots. Also, not included in the heat trans- 
fer area is the peripheral area of the slots normal to the flow, If 
the 20° skew perforated j factor were based on the total available area 
for heat transfer, the values plotted would be decreased approximately 
40%. The perforated material still represents a substantial improvement 
in the performance. This improvement is attributed to the reduction of 
inefficient stagnant areas within the matrix and the alteration of the 
thermal and hydrodynamic boundary layers caused by the slots, 

Figure 13 and Table VI present the results for a "parallel plate" mat- 
rix constructed from brass shim stock. As can be seen in Figure 13, the 
heat transfer performance falls far below that of theoretical laminar 
flow performance for parallel plates. Several things may account for this: 
the aspect ratio of the flow channels in this matrix was 77 vice ©&© for 
truly parallel plates, and also, the flow channels were not uniform in 


cross-section from one to another. The non-uniformity in cross-section 
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results in uneven mass flow through the matrix, hence inefficient use 
of the available heat transfer surfaces, and therefore a decrease in j. 

The decrease in friction factor is possibly due to the sharp angle 
regions where the formed sections contact the flat plate. The wall shear 
stress in this type of region is small which would tend to reduce the 
value of f£. These regions would also tend to produce stagnant regions 
previously mentioned, having a diminishing effect also on j. 

Figure 14 and Table VII present the results for a "parallel plate" 
matrix constructed from perforated nickel plate. This matrix was construct- 
ed in the identical fashion as the brass matrix, Although it is not read- 
ily apparent in Figure 14, the heat transfer characteristics have been 
improved nearly 80%, while the friction factor has remained the same. 
Again it is necessary to make note of the fact that the heat transfer area 
of the two matrices were calculated on different bases. These bases are 
the same as discussed for the 20° skew matrices. Applicable also is the 
approximate reduction by 40% of the plotted values of the j factor if 
all of the area, in the perforated material, available for heat transfer, 
were considered in the calculations. 

Figure 15 is a comparison of the two "parallel plate" matrices. The 
improvement realized by using perforated material is readily apparent, 
even considering the 40% reduction mentioned above. The additional per- 
formance gain in the transferral of heat is again attributed to the re- 
lieving of the stagnant areas and the aiteration of the boundary layers 
caused by the slots, 

The fact that the friction factor did not increase does not lend 
itself to previous conclusions (20° skew); that the presence of slots 


increases the friction facto:. Some additional factor must be necessary. 
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Perhaps, the skewing coupled with the slots contribute more to friction 
than straight, slotted channels. 

In the determination of the j factor for the perforated "parallel 
plate" matrix the conduction parameter, A » was considered zero, Here 
again the plotted points departed from the expected straight line indicat- 
ing that » » in fact, was not zero. Corrections for longitudinal conduc- 
tion were made using the matrix length, L, as an approximation for the 
conduction path length. These points were plotted also. Since the cor- 
rected points are overestimations and Af 0, the true value of j lies 
between these limits. 

The designer of a compact heat exchanger is generally restricted as to 
frontal area, volume, weight and pressure drop. There is presently coming 
into use, several means of comparing heat transfer surfaces with these 
restrictions in mind. Essentially, various parameters are compared 
graphically so that the relative merit of the various surfaces may be es- 
tablished. 

The ratio j/f is plotted vs. Reynold's Number as a means of relating 
frontal flow area with pressure drop. For a given pressure drop a high 
j/f£ characteristic implies a small flow frontal area. j/f when plotted 
vs. N, is known as the Flow Frontal Area "Goodness" factor. 


R 


The heat transfer power, h is plotted vs. the flow friction power, 


SED. 


orp? and the resultant characteristic curve is a measure of relative core 


volumes. Generally, the surface with a high herp vs. Earp characteristic 


will tend to have a small core volume, herp plotted in this fashion is 
known as the Volume "Goodness" factor, 
Figure 16 compares the 20° skew matrices and the "parallel plate" 


matrices. Although the "Goodness" factors are not infallible, they are 
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oF, CONCLUSIONS, 

Based upon the correlation of present data with that obtained by 
previous experimenters, the experimental test rig design presented here- 
in is considered capable of producing reliable test data on the heat trans- 
fer and friction characteristics of a core module, 

Basic heat transfer and flow friction characteristics for a new 
geometry and a new idea in material are presented in Figures 11 to 15. 

The geometric and physical properties are summarized in Figures 5 and 6. 
There appears to be a definite advantage to constructing matrices from 
perforated material, Performance is improved considerably and based on 
the ''Goodness" factors (Figure 16) size and weight advantages are realiz- 
ed, 

The 20° skew perforated nickel matrix was the top surface in perfor- 
mance, This matrix had the highest Flow Frontal Area "Goodness" factor, 
as well as the highest Volume to surface "Geodness" factor. This does not 
mean that this is the "optimum" surface for any given application but it 
does merit consideration, 

In both geometries tested (20° skew and "parallel plate") the per- 
forated material out-performed its counterpart. This would also indicate 
that further testing of this type of material (slotted) is warranted. 

It would seem that this material and some geometry should result in an out- 


standing heat exchanger surface. 
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RECOMMENDATIONS FOR FUTURE WORK. 


1) 


2) 


3) 


Some means should be found to reliably determine Le in the 

high N, range Cth Ny < 3.5). Most promising is the use 

of a periodic or cycling technique in which effectiveness is 
measured and used to find ie rather than the maximum slope. 
This ability is necessary if performance curves are to be extend- 
ed to NR = 500% 

More work can be done with perforated material. Different 
geometries should be tested to see if consistantly high per- 
formances are realized as they have been in the past. Addition- 
ally, this type of surface is suspect of fouling. If these 
materials are as good as they seem to be, then examination of 
the effects of fouling on their performance is in order. 

As more surfaces become available, they should be tested 


and their heat transfer and friction characteristics determined. 
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Matrix Material Steel 
Specific Heat (c,) BITU/1lb °F 0.12 


Thermal Conductivity (k,) BTU/hr °F ft 10.0 


Material Thickness, inches | OnOn} 
Sees feat Trenefer Aree (hmm 4.0159 
Frontal Area (A;,,) £t° 0.07055 
Total Conduction Area (A,) £t° 0.02215 
Flow Cross Section Area (AQ) £t° 0.0484 
Matrix Volume (V,,) ft 0.00199 
Matrix Density (pq) lb/ft? 155. 
Hydraulic Diameter (Dy) £t GEO0L 6 
Compactness (A) £t7/ft° 2009. 
Porosity (p) 0.686 


FIGURE 3 Physical and Geometric Properties of 
60x60 Mesh Screen Matrix 
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Matrix Material 


Specific Heat (c,) BTU/1b °F 
Thermal Conductivity (ke) BTU/hr °F ft 


Material Thickness, inches 

Total Heat Transfer Area (A) ft 
Frontal Area (Ay,.) £t° 

Total Conduction Area (Ag) £t2 
Flow Cross Section Area (A,) £+° 
Matrix Volume (V,) £t° 

Matrix Density Pa ) lb/ft? 


2 


Hydraulic oui (Dy, y ft 
Compactness (f) £t 2 
Porosity (p) 





IE 


Type 302 
DS. Steel 


Omad 
720 


0.001 
30.12 
0.07023 
0.00752 
0.06271 
0.01170 
52.9 


001388 


2573. 
0.893 


FIGURE 4 Physical and Geometric Properties of 


Plate-Fin Reference Matrix 
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Matrix Material 
Specific Heat (é, ) BIU/1b °F 
Thermal Cond. om ) BTU/hr °F ft 


Material Thickness, inches 
Total Heat Transfer Area (A) ft 
Frontal Area (A, Y £t° 

Total Riatiod | Area (A, ) et? 
Flow Cross Section Area (Ay ) £t° 
Matrix Volume ie, ) ft? 

Matrix Density Pr ) lb/ft? 


Hydraulic Diameter (D3) ft 
Compactness (6) £t7/ft? 
Porosity (p) 


2 


Type 302 


S.oteel 
0.11- 
0159) 


0.001 
18.58 
0.07024 
0.00464 
0.0656 
0.0117 
32.4 


0023581 
1570. 
0.934 


10° SKEW 20° SKEW 


202 


SKEW 


Type 302 ‘perf. 


S.Steel 
Oil L 
T.0 


0.001 
18.04 
0.0702 
0.00452 
0.0657 
0.01157 
B2e22 


002416 
1549. 
0.936 


FIGURE 5 Physical and Geometric Properties of 


Skew Matrices 
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Nickel 
0.1065 
38.7 


0.0022 
15.41 
0.0734 
0.00349 
0.06412 
0.01211 
50.1 


002528 
1272. 
0.873 











Matrix Material 


Specific Heat (e,) BIU/1b °F 


Thermal Conductivity (k,) BTU/hr °F ft 


Material Thickness, inches 

Total Heat Transfer Area (A) ft 
Prontal Area (Ap,) f£t° 

Total Conduction Area (A) rt 
Flow Cross Section Area (A,) ft 
Matrix Volume (Vv tt? 

Matrix Density (Q,) lb/et? 


2 


Hydraulic ome 7 (D,, ? ft 
Compactness (8) ft 2724) 
Porosity (p) 


2 


2 


710-30 
Brass 


0.092 
66. 


0.001 
22.057 
0.06986 
0.00551 
0.06434 
0.01164 
41.5 


~ 001944 
1895. 
921 


FIGURE 6 Physical and Geometric Properties of 


"Parallel Plate" Matrices 
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Perf. 
Nickél 
0.1065 
38.7 


0.0022 
18.586 
0.0695 
0.00416 
0.0584 
0.01158 
63.6 


001928 
1605. 
«840 
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Figure 7 60x60 Mesh Stainless Steel ‘Screens 7 
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Figure 9 10° Skew Stainless Steel Matrix 
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FIGURE 22. Matrix Holder modified with probe 
for velocity profile measurement 
and temperature distribution 
meesurement (Downstream View) 








FIGURE 23. Nichrome Wire Heater (Upstream View) 





FIGURE 2);. Nichrome Wire Heater (Downstream View) 
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TABLE I SUMMARY OF HEAT TRANSFER AND FRICTIONAL RESULTS: 
60x60 MESH STAINLESS STEEL SCREEN MATRIX 


HEAT TRANSFER RESULTS: 
en RESULTS : 
RUN REYNOLD'S NON CONDUCTION MAXIMUM N 


NUMBER ST PR PARAMETER § SLOPE vats 
1 431.57 .0267 0.0 58333 3510 
2 234.81 60425 0.0 668566 4.96 
3 159.26 20521 0.0 -74639 6.09 
4 129.69 .0620 0.0 ~80395 725 
5 106.78 .0659 0.0 ~82926 Ton 
6 79.88 .0755 0.0 -87817 8.83 
7 50.69 .1040 0.0 1.05579 11.90 
8 29.60 01535 ise 1.22805, 17.94 
9 15.43 ©2584 O-G) 1.57029 30.19 
10 — 8.49 -4760 0.0 2.11897 55.65 


ISOTHERMAL FRICTION RESULTS: 
ISOTHERMAL FRICTION RESULTS 
RUN REYNOLD'S FRICTION 


# NUMBER ‘FACTOR 
i 437.93 -68548 
2 238.27 ~82196 
3 161.60 ~89263 
4. Dsies ©97400 
5 108.36 1.06410 
6 81.07 1.23243 
7 51.44 1.57979 
8 30.04 2.31547 
9 15.66 4.47090 
10 8.62 9.50629 


Sy! 





TABLE II SUMMARY OF HEAT TRANSFER AND FRICTIONAL RESULTS: 


PLATE-FIN REFERENCE MATRIX 


HEAT TRANSFER RESULTS: : 


RUN REYNOLD'S NAN 


# 


CWO OA AU FU DH 


~ 


NUMBER 
Np 
505.17 
453.26 
387.19 
322.08 
248.92 
187.25 
125..74 
75.04 
42.68 
23.81 


siNpae/? CONDUCTION MAXIMUM 
PARANETER SLOPE 

J 
00277 =» .00128 «52583 
00429. = 00150 «55712 
00505 00167 «=. 58044 
00634 .00201- 62383 
00842 ~00260 693577 
.01162 00346 19434 
01761 .00515. «94678 
02889 = .00862,-—s« 1.16341 
~05699 ~01516 1.43365 
08403 02717 «1.43851 


ISOTHERMAL FRICTION RESULTS : 


RUN REYNOLD'S FRICTION 
FACTOR 


# NUMBER 
512.59 
439.62 
392.86 
326.79 
252.57 
189.99 
127.58 
76.14 
43.30 
24.17 


wo OA AW FUND 


j~ 
o) 


f 


02197 
02485 


«02795: 


.03304 
404210 
»05433 
.08205 
~12974 
»23164 


039637 
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1.GS: 


2.60 
3.06 


3.84 ° 


5.10 
7.04 
10.67 
18/ Gal 
34.54 
50.93 


(A=0) 


1.68: 
2.60 
3.06 
3.84 
5.10 
7.04 
10.44 
16.21 
25.07 
25.22 


TABLE III SUMMARY OF HEAT TRANSFER AND FRICTIONAL RESULTS: 
10° SKEW STAINLESS STEEL MATRIX 


HEAT TRANSFER RESULTS: | 


RUN REYNOLD'S N..N_.°/? 


# 


wow On nn fF UP FH 


NUMBER 
Np 


805.27 


643.32 


517.86 © 


402.58 
305.00 


200.72 


68.74 
3761 


j a 
.00375_ 00081 
00440 00101 
00695 00126 
.01073 .00162 
01488 00214 
02271 00324 
03642 200544 1 
.06376 .00953 1 
.10788 01728 


MAXIMUM 
SLOPE 


© 46533 
.51088 
255129 
62114 
- 70167 
83894 
.02840 
-27448 
44693 


ISOTHERMAL FRICTION FACTOR: 


RUN REYNOLD'S FRICTION 
FACTOR 


# NUMBER 
“ae 
1 817.04 
2 652.71 
3-25.43 
4 408.46 
5 309.46 
6 203.79 
7 121.41 
8 70.00 
9 38.16 


au 


~O2T91 
203405 
04092 
004951 


06345 


209412 
015138 
«27472 
47276 


1.32 
os 
2.45 
3.18 
5224 
8.00 
12.83 
22.46 
38.00 


tu 

(A =0) 
Ley 
i255 
2.45, 
3.78 
5.24: 
7.98 
L246 
19.62 
25.59 





TABLE IV SUMMARY OF HEAT TRANSFER AND FRICTIONAL RESULTS: 


20° SKEW STAINLESS STEEL MATRIX 


HEAT TRANSFER RESULTS: . 


RUN REYNOLD'S NomNpp 


att 


OWN MYnAWU UD 


NUMBER 
a 

835.38 
728.65 
636.57 
535.10 
412.85 
312.06 
208.65 
126.47 
71.21 
39.79 


00394 
00414 
00475 
.00805 
01182 
01524 
02167 
03676 
~06593 


09441 


2/3 GonpUCTION MAXIMUM 
PARAMETER § SLOPE 


A 
.00078 47359 
00090 48768 
.00103 52394 
.00122 56474. 
.00158 63657 
.00209 70256 
.00313 81581 
.00517. =: 1.02339 
.00917 1.28940 
01646 1.39790 


ISOTHERMAL FRICTION FACTOR: 


RUN REYNOLD'S FRICTION. 


# 


i 
2 
3 
4 
5 
6 
7 
8 
9 


10 


NUMBER 
NR 
848.11 


139.39 
645.95 


542.97 
418,94 
316.65 
211.73 
128.33 

72.26 

40.40 


FACTOR 
ag 
03197 
005471 
-03908 
04465 
005431 
06856 
10255 
16015 
+ 28784 
. «48428 


We eis 
1.43 
1.64 
2.78 
4.09 
5.26 
7.48 
12.69 
22.76 
32.59 


(A=0) 


1.56 
1.43 
1.64 
2.10 
4.09 
5.26 
7.48 
Oe 2) 
20.09 
25.79 





TABLE V SUMMARY OF HEAT TRANSFER AND PRICTIONAL RESULTS; 
20° SKEW NICKEL MATRIX (PERFORATED PLATES ) 


1 
2 
3 
4 
5 
6 
7. 
8 
9 


10 


HEAT TRANSFER RESULTS: 
EA RESULTS 3» 


RUN REYNOLD'S NooN,.°/> conpucrron 
# NUMBER . - PARAMETER 
NR j d 
966.07 00900 00337 
835.43 .01146 -00390 
738.42 -01219 00441 
619.79 .01418 200526 
484.37 .01802 00674 
362.97 .02497 .00898 
239.79 203517 201359 
145.03 205199 202243] 
82.26 07879 03951 1 
46.22 .08147 -07039 1 
23.50 .09134 -13843 od. 


‘121 


MAXIMUM 
SLOPE 


~ 56394 
-60157 


e 61776 BS 


- 64762 
71212 
81813 


095653 


~ 14669 
- 29808 
42079 
+ 50225 


ISOTHERMAL FRICTION FACTOR: 
a LUN PACTOR 


RUN REYNOLD'S FRICTION 


NUMBER 
NR 

1 980.29 

2 847.72 

3 749.28 

4 628.91 

5 491.63 

6 368.31 

T 243,32 

8 147.16 

9 83.46 

10 46.89 
11 23.85 


FACTOR 
re 

- 03316 

~03709 

-04116 

-04686 


-05802 


.07208 
10258 
16737 
29667 
49889 

1.09090 
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Zane 
3.42 


Bde) . 


4.24 
5.41 
7.58 
11.29 
19.65 


tu : 
(As0) 
2 ie 
3.46 
3.68. 
4.28 
5.44 
7.54 
10.62 
15.70 
23 ao 
24.60 
27.58 





TABLE VI SUMMARY OF HEAT TRANSFER AND FRICTIONAL RESULTS: 


"PARALLEL PLATE" BRASS MATRIX 


HEAT TRANSFER RESULTS: . 


SBE 


MAXIMUM 
SLOPE 


°49410 


54861 | 


259598 
«66717 
015472 
~87166 
~97562 


RUN REYNOLD'S mene °/ 0 CONDUCTION 
NUMBER PARAME 
Np j 
683.97 .00376 00896 
551.23 600527. ~—.01112 
437.55 00748 01402 
340.15 01038 .01802 
255.77 01335 .02390 
168.67 02148 03634 
202.16 03483 05999 
57.06 = .10761 1.07854 
= 18369 


98829 


ISOTHERMAL FRICTION FACTOR: 


RUN REYNOLD'S PRICTION 
FACTOR 


# NUMBER 
NR 
1 694.19 
2 559.46 
3 444.09 
4 345.24 
5 259.52 
6 171.19 
7 e327 
8 57.90 
9 33.86 


a 


02846 
005454 


04292 


05240 
~06907 
- 10680 
~17239 
032370 
© 53085 


N N 


tu. Gia 

(A=0) 

1.44 1.46 
2.27 2.40 
5.22 Biwi 
4.47 4.63 
5 Bee) 
9325. Gav 


15.00 11.15 





TABLE VII SUMMARY OF HEAT TRANSFER AND FRICTIONAL RESULTS: 


"PARALLEL PLATE" NICKEL MATRIX (PERFORATED PLATES) 


HEAT TRANSFER RESULTS :. 


RUN REYNOLD'S NN 


# 


oO OA AWM FY PY FH 


NUMBER 
NR 
819.49 
647621! 
518.89 
400.50 


304.49. 
200.68 


121.15 
67.53 
3792 


2 
J d 

00590 «00393 
«00810 00497 
00988 00620 
01458 00804 
01891 01058 
°02659 01602 
~04039 -02658 
205502 -04768 
004945 


08478 


MAXIMUM 
SLOPE 


54 (em 
»58864 
«62923 
eioleo 
081941 
°95459 
1.16142 
ipvo4a SO 
Pe2ngol 


ISOTHERMAL FRICTION FACTOR: 


RUN REYNOLD'S FRICTION 
FACTOR 


# 


1 
2 
o) 
4 
5 
6 
7 
8 
9 


NUMBER 
NR 

831.56 
656.72 
526.51 
406.49 
309.04 
203.63 
122.96 

68.54 

38.48 


ai 
02395 


202955 


.03577 
-04439 
.05682 
.08725 
»14242 
»26776 
45874 
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2.52 
9.\a 
35.89 
5.77 
7.61 
1 9k 
21.61 


tu 
(A=0) 
2.56 
3.24 
Bosh) 
veh). 
7256 
LOgoo 
1695 
22.00 
19.70 
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APPENDIX A, 
DATA REDUCTION RELATIONS 

The following summary of some of the data reduction relations is 
given to clarify some of the definitions of the geometrical parameters, 
to show how they influence the deduced results and to facilitate an 
understand as to the exact procedure required to produce from raw data 
a point on the j vs N, and f£ vs Np curves, 

The three most significant geometric parameters that afford a means 
of comparison between matrices are so defined that there exists a unique 
relation between them. Independent measurement of any two establishes 
the third or serves as a check on its independent measurement, The 
three primary surface geometry parameters of a matrix are: 

1) HYDRAULIC DIAMETER: 
ac, = 4G 1 (A-1) 


2) POROSITY: Matrix Flow 
Void Volume ALL A. 
p = Total Volume = —*—_ = -<* (A-2) 
Apes A 
f-% fr 


3) AREA COMPACTNESS: Total Heat Transfer 
Surface Area 
p = Total Volume = (A-3) 





NOTE: If equation (A-2) is divided by (A-3) the unique relation mentioned 


previously is formed. 


a P/B (A-4) 


MAXIMUM SLOPE: 


As previously discussed, through the work of Locke, Mondt, Crestwick 
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and Howard, it has been established that the maxi slope of a general- 


ized heating or cooling curve is a unique function of Ne and A , the 
ui 


conduction parameter, 





peat 
OT Crammer e 
a 4 = £{(N oh 5 
d(7/N.) ) | max a 
tu 
where 
‘T = generalized time variable 
h A =m t¢yi#tial 2 
T = eC t, = initial matrix temp. 
ss 
ral = air temp. upstream 
and N _ ge t = air tem downstrea 
ou fe po = at p. do am 


therefore: 


= ; “f 
aN ne (We, © 


and 4( T/N.) = me UPTT GlS)9 Sponeepeosnnnenaicacenncyid: 3005 (25) 





also: 
ed 7 See a) Deo 
al” Ps Cet at 
t e 
fe i 1 
and: d¢ i ) = = df t = t ye eoeeove (A-6) 
try t. ey ct, £2 fal 
combining (A-5) and (A-6) 
t -t¢t 
£2 al 
As E Woes 1 SEG) TSA 
oe i 2 ao 7) 
d¢ a7 Ne max P faa i 


The matrix is heated to a uniform temperature, tie The heat source 


is secured and the difference between downstream and upstream temperatures 
(t eo - try) is recorded vs. time. The maximum rate of change of this 
difference with time may be determined from the record. 


From Figure 17: 


d(teo - Coy) a 
d9 max 


De 
x 
x 
Chart speed seccar sce 
y= d(t ro - te) inches 
MS (tery = t,) inches 
WoC, Matrix heat Capacitance _ CT. BTU/°F 
a Fluid heat Capacitance — Aas BTU/°F sec 
Pp 


Finally, Equation (A-7) becomes; 


Ag Se) ES 
Ss 


dC FTN, mie . . Chart speed 


With this value of the maximum slope and a calculated value of A ; 


Howard's curves or tables may be entered to obtain a value of Neu 
MASS FLOW RATE: 
The mass flow rate is calculated as outlined in reference (7): 


m= 359 KD? Fay [ar, % 
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where K = flow coefficient, velocity of approach factor 
Fa = thermal expansion factor 
Y = expansion factor 


D = orifice diameter, inches 


° 


y = specific weight of fluid flowing at the orifice 
temperature ye corrected to existing 
atmospheric pressure 

AP, = Pressure drop across theorifice , inches Ho 
REYNOLD'S NUMBER: 


By definition, the Reynolds number is 


N, = - t— (A-8) 


Where G is the mass flow velocity based on the flow mean cross- 


sectional area, Aa: 


G= a Geooa = PT Gene " (A-9) 
c P “fr 


Note, if (A-9) is substituted into (A-8) 


N. = 


ae aa A-10 
2° ae Cape (A-10) 


; y, gilli 4 i 4m il 
but from (A-4) H/p ao N= A. (2) (A-11) 


= B . v7 
Therefore, for a given range of mass flows, the Reynolds number range 


for a given matrix is proportional to the inverse of its compactness, 


No 1B 


For ease of calculation, recall from (A-3) B = + 
fr 


substituting into (A-11) yields 
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THE FANNING FRICTION FACTOR: 


The following sketch describes the flow system to be considered. 





al a b 2 
VILLLLLILLLILLLIL LL LL LL 
ee ee , ies 
ane CLLLLLLLLLILLLLULLLLD a 


LLLLLLLLLLL ALLEL TLL LL. 


For gas flow heat exchanger application, the pressure changes 
from section 1 to a and from b to 2 are very small, relative to the 
total pressure, and therefore Vax Vy and VU, = Va 

The relation for the flow stream pressure drop calculation for 


most heat exchanger cores is (8): 


P, 29 P, 


S22. G vi Chest pt)+ 2(H -1) +f A Yes = Cage ees (A-13) 


entrance flow core exit 
effect acceleration friction effect 


Solving (A-13) for £ and recalling from (A-1) A/A, = L/r,, also 


Ve= e » we find: 
aly oy | Mh 
fz |, (8) -( ta Pm ~ fom (A> (I+ PY ‘i (A-14) 


Now the perfect gas law may be applied since testing is done with air at 
moderate temperatures and pressures. Equation (A-14) becomes for the 


isothermal case: 
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5 
i€ ' ‘ e 
} 
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> = —_ = 
4 « = a 
i; 6 
7 => 


= [2 QE) EEK) EAP) | ces 


where e = Pp ane 9 er 


The first term is by far the largest contributor to the friction factor 


for small pressure differentials, therefore making the approximations 


Equation(A-15) becomes 


f= | Pol) Gk de aB (1p) gets 


where: Ke entrance coefficient 
Ke = exit coefficient 
whose values are determined from appropriate curves in Chapter IV, refer- 
ence (8). 
This is the equation used to determine the friction factor as 


presented in this thesis. 


If a further assumption is made such that 


= 240 Pn SB (A-17) 
f Se frm at sa 
and equations (A-4) and (A-9) are substituted, Equation (A-17) becomes; 


f fa PeOP (Ae) (A-18) 


rh 4 
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Therefore, for a given matrix the friction factor is directly pro- 
portional to the cube of the porosity and indirectly proportional to the 
compactness, 

£ > P/Q (A-19) 
COLBURN j FACTOR: 

We have A and maximum slope with which we get Ne from Table VIII 

or Figure 18 taken from reference (5). 


The j factor is defined as: 








h 
j = Non Nop 2/3 e x Ae 273 (A-20) 
substituting (A-9) for G and multiplying by Ts 
hA Ac 2/3 
j= fe A pr Cane 
by definition: J hA 
tu mCp 
Ac Zi 3 
Therefore: j= Ne A Nor (A-22) 
Note, from (A-1) Ac/A = "a/L, substituting P/2 for 
th from (A-4) gives: 
, 2/3 1 
j=™ Noe . P18 (A-23) 


which leads to the relation that for a given matrix, j is directly pro- 


portional to the porosity and inversely proportional to the compactness, 
i oc %e 
HEAT TRANSFER POWER AND FLOW FRICTION POWER: 
The heat transfer power per unit area per degree temperature differ- 


ence is derived as follows: 
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J = Noor 
ox gp eee K 1 
a i N 
ST NRNPR k - yz R 
i sian 1 ee 
es NR PR 
Solving for h: 
“5 Md 1 
hee =e Se eg N. j 
‘ 2/3 Diy R 
IPN 


p are evaluated at standard properties which by 


convention are taken as those for dry air at 500°F and one atmosphere; 


Now if cx and Np 


c 0.2477 BIU/1b °F 


Pp 
AA = 0.0678 1b/hr ft 
fO = 0.0413 1b/£t° 
Nop = 0-671 


Thus, the heat transfer power per unit area per degree temperature 


difference evaluated at the above standard conditions, Dorp? becomes 





he OL Oo (N, 5) BTU/hr £t2°F 


STD Ar, 


The flow friction power per unit area is derived as follows: 


aes panels (cay 
[ 


£ cz 


AP( 2) = ro 
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E =ap(S (aa) © Fit 5. (4 mea) 


where Se rP V 


> #1 a) (Bi) Ne 


2% 


Now if and f are evaluated at standard properties (previously 
defined) the flow friction power per unit area evaluated at the above 


d i : 
standard conditions, E orp? becomes 


N 
a -7 i 3 R 3 [ a 
Earp = loll 10 ar, (500) ££ \hp/fe 


Figure 16 provides a comparison of the 20° skew and the "parallel 
plate" matrices from a heat transfer power vs. friction power point of 
view (hoon vs. orp)’ For comparison purposes the surface geometries 

=o 


have been adjusted to a common hydraulic diameter of At, = 2x LO MMEE. 
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APPENDIX B. 
DESCRIPTION OF EQUIPMENT 

To obtain heat transfer data, a transient technique was employed. 
Briefly, this method consists of heating the test matrix, to a uniform 
temperature and then subjecting the matrix to a step change in fluid 
flow temperature. The fluid temperature downstream of the matrix is 
recorded versus time. 

Friction factor was calculated from pressure drop data obtained 
from static pressure taps located in the test section immediately up- 
stream and downstream of the matrix, 

The equipment necessary to complete the experiment falls into one of 
the following systems. 

1) Fluid source 

2) Flow metering system 

3) Temperature measuring system 

4) Pressure measuring system 

5) Fluid heater 

6) Matrix holder & test section casing 
FLUID SCURCE: 

The working fluid was, of course, air, Air was provided to the 
apparatus by placing the entire rig on the inlet to a 100cfm, 6 stage 
centrifugal air compressor. The compressor requires an operating voltage 
of 120” at 90amp. 

FLOW METERING SYSTEM: (See Figure 1) 

The flow metering system consists of an orifice in a 3.08" inside 

diameter aluminum tube. The orifice section was constructed in accord- 


ance with the ASME specifications of reference (7) for standard D and 


D/2 pressure taps. A large flow range was obtained by means of 7 orifice 
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plates of throat diameter 2,310, 1.971, 1.232, 1.540, .462 and .308 
inches, respectively. Accurate control of the pressure drop across the 
orifice was maintained by a large sliding valve downstream of the ap- 
paratus on the inlet to the compressor and another sliding valve on the 
discharge side of the compressor. There is an additional valve on the 
inlet side of the compressor that is fully open or fully closed. It is 
normally opened at flow rates less than 150 lb/hr to prevent stalling the 
compressor. 

TEMPERATURE MEASURING SYSTEM: (See Figure 1) 

Without exception, all temperatures were measured with iron-constan- 
tan thermocouples. A single thermocouple measured the orifice tempera- 
ture a) just upstream of the orifice at the inlet to the apparatus. 
This temperature was read in millivolts on a RUBICON CO. portable pre- 
cision potentiometer, converted to temperature and recorded for each run. 
Measuring Ty is a group of 5 thermocouples in series, bound together and 
insulated from each other by teflon tape and inserted in a 1/8" diameter 
aluminum tube, The tube was placed in the air stream with a small window 
cut in it facing upstream, At this point ambient air temperature was 


measured. T, and T, are 5 wire thermocouple grids permanently mounted in 


2 4 
the test section cabinet. Their wires are connected in series to magnify 
the thermocouple output voltage. T, is a movable 5 wire thermocouple 
grid located in the matrix holder and is placed adjacent to and down- 
stream of the matrix, it's wires are also in series. It's position is 
apparent in Figure 21. The thermocouples arranged in this manner provided 
flexibility and complete coverage of desired points of interest. If the 


thermocouple output voltages from Ty and T, are bucked against each other, 


ih 
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the temperature measuring system response curves may be generated on a 

Leeds & Northrup 'Speedomax" temperature recorder, T, vs Ty, gives the 
indication of temperature uniformity across the matrix and T, vs T, 

records the rate of change of the difference between upstream and down- 
stream temperatures with time. It is this record that provides the pri- 
mary data for the evaluation of heat transfer data. The precise relation 

of this record to the Colburn j factor is explained in Appendix A. The 
record is made on a Leeds and Northrup “Speedomax" Temperature Recorder. 

The recorder has variable chart speeds, which permits maintenance of accept- 
able cooling curves with regard to slope determination. The Speedomax also 
has a range of 0-20 mv. For this experiment, the recorder is precalibrated 
on a 3 mv. scale. 

There is one addition thermocouple grid that is not used in the experi- 
ment, but was used in checking out the apparatus. This grid would take the 
place of T, and its purpose was to determine the temperature distribution 
across the cross section in both the vertical and horizontal directions. 

It was placed in a special holder that also held a velocity probe (See 
Figure 22) and was lined with balsa wood to provide a smooth and continu- 
ous internal flow passage. 

PRESSURE MEASURING SYSTEM: (See Figure 1) 

Pressure taps are located upstream and downstream of the orifice and 
upstream and downstream of the matrix. Each tap is connected by Saran 
tubing to an appropriate manometer or draft gage. The following instru- 
ments were used: 

1) Type "C'' micromanometer, E. Vernon Hill & Co. 0-1,25' 
2) Ellison Differential Direct Draft Gage, 0-6", 0-8" 


3) Ellison Inclined Draft Gage, 0-6" 
4) Meriam Instrument Co. 0-60", 
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Any one or combination of these instruments can be used to measure 
the differential pressures of the orifice and the matrix. Frequent cross 
checks were made to ensure reliable operation. 

FLUID HEATER SYSTEM: (See Figure 1) 

The air heater section consists of 28 nichrome .0031" diameter wire 
heaters. The heater system is designed to heat 1000 pounds of air per 
hour, 20°F above ambient temperature. Nichrome was chosen because of 
its high resistivity, low thermal conductivity and low specific heat, 

The nichrome wire is wound on a bakelite frame, two heaters to a 
frame. The wires are 1/32" apart and there are 50 and 52 wires respective- 
ly, to each heater on a frame. The heaters on each frame are connected in 
parallel and then connected through a switch, to a variable voltage bus. 
There are two frames to a plastic holder, as shown in Figures 23 and 24, 
In these figures, the heater frame and its electrical connections can be 
clearly seen. 

A schematic wiring diagram for one frame is shown below. All 14 
frames are wired in this manner, putting all 28 heaters in parallel with 


the capability of switching on or off two at a time, 





Bus BAR 

ee 
SELECTIO 
PANEL 


a 


oFF- ON 
Sswitca 


SELECTOR 
SwitcH 
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The 28 heaters arranged in this way, provide a variety of flexi- 
bility. The number of heaters in use, should decrease as the flow rate 
decreases, Voltage variation can be obtained by a GENERAL RADIO COMPANY 
"VARTAC", 0-240”, 8a, 50-60 cycles. It is recommended that a voltage of 
220” or in this vicinity, be maintained, at all flow rates, and the 
number of heaters used, decreased until the number of heaters is four. 

At this point, the voltage should be decreased, as necessary, to obtain 

a 20° temperature rise, the number of heaters remaining constant. A 
note of caution: Do not energize a heater, unless the air is flowing, 
otherwise, the wires will sag due to thermal expansion and short circuit- 
ing will occur. 

MATRIX HOLDER AND TEST SECTION: 

The matrix holder and test section casing are made from polyethylene 
plastic. This material was chosen because of its low heat capacity. 


The test section casing holds thermocouples T, and T, permanently 


2 
mounted upstream and downstream respectively. It also contains the up- 
stream and downstream static pressure taps used in determining the fric- 
tion factor. The parts making up the casing were machined to close toler- 
ances, to ensure a tight fitting matrix holder, and a minimum of air leaks. 
The parts comprising the matrix holder, as implied above, were machin- 
ed to close tolerances to ensure a tight fit. Great care was taken to 
ensure smooth alighment of the flow channel through the heaters and the 
test section. The size of the matrix holder was arbitrarily chosen. 
The flow channel dictated the width, and the length is such that matrices 
of various flow lengths, up to 3", may be tested. The flow channel is 


3-3/16" x 3-3/16". 


The matrix holder has one special provision and that is for the movable 
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thermocouple T A small hole is drilled in the front, and the wires 


ee 
are led through. The matrices are placed in the holder and held in 
place by bits of Styrofoam plastic cut to the size necessary to hold 
the matrix square in the flow channel. 

That completes a description of the various systems involved in 
the apparatus, but there is one other item that bears mentioning. Just 
upstream of the heater section, is a section provided for screens to 
straighten the velocity profile. Provisions were made for 1 to 4 screens. 


One 60 mesh stainless steel screen was found to be sufficient, to straight- 


en the flow at the flow rates used (40-1000 lb/hr.) 
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APPENDIX C 


CALCULATION OF GEOMETRIC CONSTANTS 
FOR PERFORATED MATERTAL 


As demonstrated in Appendix A, the surface geometry has a strong 
influence on the heat transfer and friction factor curves, Specifically, 
the effects of p and B were examined. 

In using perforated materials, the presence of the slots, require 
modification of these two parameters. It is the intent of this Appendix, 
to demonstrate the calculation procedure utilized to determine the geo- 
metric parameters, on which the graphical and tabular results for the 
matrices constructed from perforated material are based. 

The material used in the perforated plate matrices, was a pure nickel, 
type 160/40TV, slotted opening plate material, .0022 inches thick. This 
material is a product of Perforated Products, Inc. 

The basic geometry of the plate is: 


+ O22¢m 


fo -04zem — ie ke :] 
le .o38con ale . o2dons i 





Using the average between the top and bottom of the tapered slots, 


consider an average slot pattern as: 


— 023m <—_.p40em —> 
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CROSS HATCHED AREA Le = (.063) (.016) = .001008 ene 


AREA OF SLOT IN Aj, = (.040)(.0020) = .00008 oe 


Now define ‘i = Aslot = AREA REDUCTION DUE TO SLOTS 
A 
ch 
90008 } 
Y ~  ,001008 — ad 

Therefore, plate Porosity (p,) = (1 - Yo») = .9207 
Meu define S$ _ CROSS SECTION AREA OF SLOTS _ CONDUCTION AREA 

ee CROSS SECTION AREA REDUCTION DUE TO 

SLOTS 


(040) (.0022) _ ; 
& = (.063)(.0022) = 03+5% 


The following calculations are for the 20° skew perforated plate 
matrix: 


DIMENSIONS OF PLATE: 4.285" x 1.98" x .0022" 


NUMBER OF PLATES: 142 
DIMENSIONS OF MATRIX: 3.255" x 3.248" x 1.98" 
WEIGHT OF MATRIX: 275.5735 grams 


MATERIAL CONSTANTS; fa = 555 1b/-,3 


k= 38.7 BIU/, op ft 


cot 1065 BTU/ 1b°F 


PLANE SURFACE AREA = 142 “ss 1.98) (2 = 16.733 £t2 


TOTAL HEAT TRANSFER AREA (A) PLANE SURFACE AREA - AREA OF SLOTS 


(1-3) SURFACE AREA 


(.9207) (16.733) 


A = 15.41 fer 


FRONTAL AREA (A,_) = 3.255 x 3.248 = 10,572 ne 


A = 0/3 et? 
fr 
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FREE FLOW AREA = FRONTAL AREA - CROSS SECTIONAL AREA OCF METAL 


= 10.572 ane - 4.285 x .0022 x 142 = 9,233 oe 


| A. = .06416 et? | 


CONDUCTION AREA (A.) (142) (4.285) (.0022)(1- & ) 


(142) (4.285) (.0022) (.375) = .50198 in” 


It 


A, = .003486 a 


MAXTIX VOLUME (V_) = 3.255" x 3.248" x 1,98" = 20,933 in? 


MATERIAL VOLUME (V,) = 142x4.285"x1.98"x.0022" = 2,650 fae 


WEIGHT OF MATERIAL (We) = 275.5735 ems./453.6 em/lb = .6075 1b. 
3 


MATERIAL VOLUME (T_) = Bie = ,6075 1b/555 lb/ft? = .001095 ft 


The porosity of a matrix has been defined as: 








_ MATRIX FLOW VOID VOLUME _ 7 Re). Vs 
7 MATRIX VOLUME ~ om s Te 
Vm 
001534 
P re 1 ra 01211 = 8733 
ve 001095 
p=l- 7a il = ‘ODL & 1 - ,09042 = .90958 


The compactness of a matrix has been defined as 


_ TOTAL HEAT TRANSFER AREA 
8 z MATRIY VOLUME 


Here, we will redefine le to include the effect of area reduction: 


_ (PLANE SURFACE AREA)(1- % ) 
eC y MATRIX VOLUME 
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VAS 
~ O67 330 0pes fte/ ft 
= O12 11 = 1381.75 


For Y= .0793 
se 2 
_ (16.733) (.9207) _ ft (ee 
6 = Sass 1272.18 ft 


The hydraulic diameter is conventionally defined as: 


a i FLOW CROSS SECTION AREA 
H WETTED PERIMETER 


Appendix A shows that hydraulic diameter can also be expressed as: 


ewee! 
D, = 4 6 
To be consistent with the conventional definition, for the case of 


a perforated material, p is used vice B » since the slots have no 


effect on the hydraulic diameter, 


Ses | | 
Dy 1381.75 .002528 ft 


Finally then, 


it 
> 
i) 
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APPENDIX D 
SAMPLE CALCULATIONS 
The following sample calculations are based on recorded data from 
Run #3 for 20° skew matrix constructed from perforated nickel sheets. 
This is chosen primarily because the basic seometric parameters are com- 
puted in Appendix C. 


RECORDED DATA; 


be = 2.310% P = 29,963 in. Hg 
[0] atm 
Are = 2.61 in. H,0 CHART SPEED = 4.0 in/sec 
AP = 3.96 in. HO W 2.6075 lb 
m z m 
P = 4,69 in. HO 
s 2 
T = 68.0 °F 


0 
MATERIAL CONSTANTS: 
co = .1065 BTU/1b °F 
ke = 38.7 BTU/hr ft°F 


PRECALCULATED GEOMETRIC PARAMETERS: (Re: Appendix C) 


A = 15.41 ne f =) ie oks 
2 A ZS 
Ac, = -0734 ft ie) = 12720 18 9te ct 
A. = .06416 ae ty = .000632 ft 
B, = .003486 £t° 
L = 1.98 inches ( — denotes corrected for slot effects) 


DETERMINATION OF MASS FLOW: as outlined in Reference (7) 


E Jarer. 
m = 359 K Do Ee né A Po ¥ 


1.0 Fig. 40b Ref. (7) 
1.0 Fig. 38 Ref. (7) 


where: Ye 


rv 


84 





. 





(359) (2.31) kK \ (2.61) ( seo) (.07522) 


m = 
m = 849.40 K 
assume K = .749 p= d/D = oa =o 
Then wh = (849.40)(.749) = 636.21 lb/hr 
Reynold's number for the pipe ND = ae 
N te (4 636.21)(12 = 7.35 x 10° 


Rp (7) (3.08) (.04410) 


From Table 5, Ref. (7) P= Obie SR ec 1, ee ee 


Therefore 5 = 636.21 Ib 


/hr 
MATRIX REYNOLD'S NUMBERS 

There are actually two Reynold's numbers to be calculated for the 
matrix. One is for isothermal flow and is associated with the friction 
factor, the other is associated with heat transfer. They differ only in 
the temperature at which the absolute viscosity JA is evaluated, Xt « 


(friction) is evaluated at ambient temperature (T.)- KX (heat transfer) 


is evaluated at ie + 10°F (an estimate of the average bulk fluid tempera- 





ture). 
4um (4) A228) £¢(636.21)? hr 
Ne * ial = 12 Ey oiseee 
£ (15.41£t2) (.04405 1b/ft hr) 
_ 4im Ae _ 04405 


CONDUCTION PARAMETER 
ne eae _ (38.7) (.003486) (12) ae 
" Tom CG)  (1,98)(636.21)(.24) * 
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MAXIMUM SLOPE 
From Figure 26 the slore of the cooling curve was determined as 
described in Appendix A. The slopes of three curves were determined and 


the average of the three was used in the calculations. 


al c 
mS 20 _ -&§ , SLOPE xCHART SPEED 
C dO / max Cc 
C, = Wc, = (.6075 1b)(.1065 *7/1b°F) = .06470°"/°F 
re . 1b BIU /5 Jos eee ee 
C=mec = (636.21 “/hr)(.24°° /1b°F) E59 -Gec 
= 04241 BT /sec oF 


SLOPE = J l6eiau ines \(ceeinteure 26) 


Finally: 


_ (.06470) (.10775) 4) - [.65745. | 
MAXIMUM SLOPE = (“575415 = | .65745 


COLBURN j FACTOR 





With x and MAXIMUM SLOPE enter Table VIII and find Neu = 4,40. 
Linear extrapolation will give desired accuracy. 


Now, A 2/3 


y Non me nee from Appendix A, Eq. (A-22) 


Bice Forel, = ad (.796) (4.40) = [ .01458 | 


FANNING FRICTION FACTOR 


As derived in Appendix A Equation (A-16) the friction factor is: 
= eB: AP\ _ = — AP 1 ta 
re pol) Gs) ale (+p) |B 


pee = seeomiagHyO x 5.204 Fo oo oie peme ee I Sipel 


sf 
2 in H,0 
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AP, = 4.69 in H,0 x 5.204 Bf Senate e695) pai 


DiS Lee 
fe in H,0 


2 
Pp. = 29,963 in H. x 040020 ta Seed 
atm g in He 
Po =P. = P = 14.7178 seeleos eeu leteesenet 
1 atm s 


\ ee oa Pe = 14,5483 - .1431 = 14.4052 psi 











2 1 
P.+P 
ee Soe 
Pele = 5 = 14.4768 psi 
_ Pmean? MM) (14.477) (144) ogyo 1p ee3 
(mean > rer. = Geeortsze) 
m _ 636.21 tb/hr _ Qo 1b 2 
G = ee Saar PNG) by los 0%: = 2.754 /sec ft 
. 000632 
Bit, | egaya  ee 
K, = -.738 K. = .949 from Chapter IV Reference (8) 
2 2 
(l+p)= 1 + (.8733) = 1.7626 
AP 
2 m _ (64.4) (.07407)(20.61)= 12.9579 
@ mean ¢2 7.587 
AP 2,  .1430eped _ 
: (l+p") = Ways psi (1+ 761) = 01742 
mean 


£ = (12.9579 - .01740 - .211)(.00383) = (12.7295) (.00383) 


-| 04859 | 


All thermal properties CL: Nop? Kg? and ce) vere obtained from 


References (12) and (13). The temperature at which the properties were 
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APPENDIX E 


A FORTRAN COMPUTER PROGRAM 
FOR DATA REDUCTION 


A short computer program was devised, Figure 27, so that from basic 
experimental data, the end result is: 
1) Conduction parameter (A) 
2) Maximum slope 
3) Reynold's number for heat transfer 
4) Reynold's number for friction 
5) Fanning friction factor 

It was anticipated that this program along with a subroutine would 
ultimately produce Ne and subsequently the Colburn j factor. 

Time did not permit the adaption of Moreland's work (11) to a sub- 
routine but the program is still useful. It's value when data for a 
single matrix is to be reduced is nebulous. However, the reduction of 
data for multiple matrices can be accomplished with a considerable savings 
of time. 


The following is a list of inputs necessary in Fortran and English 


language along with the required dimensions: 


CARD NO. 

30 XML Matrix length ft 
AS Solid conduction area fee 
RH Hydraulic radius Et 
POR Porosity - 
AFL Flow cross section area aes 
AHT Matrix heat transfer area ft? 
WM Matrix weight lbs 
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CARD NOC. 


40 EK 


CK 


SK 


CM 


N 


10 DO 


DELPO 


GAM 


ATMP 


COR 


UFR 


10+2 CS 


NELH 


HS 


SLO 


TEMPO 


The output of this program in Fortran and English 


Exit flow coefficient 

Entrance flow coefficient 

Solid thermal conductivity 

Matrix material specific heat 
Number of runs 

Diameter of orifice plate opening 
Pressure drop across orifice 
Specific weight of fluid 
Atmospheric pressure 

Orifice correction K 

Fluid viscosity for heat transfer 
Fluid viscosity for friction 
Chart speed 

Pressure drop across matrix 


Static pressure matrix inlet 


Slope of generalized cooling curve 


Orifice temperature 


their dimensions is: 


EMDOT 


RNUMP 


RNUMH 


CAPFL 


CAP‘ 


CONDPAR 


Mass flow rate (tm) 


Reynold's number for the pipe 


Reynold's number for heat transfer 
Thermal capacitance of the fluid (C) 


Thermal capacitance of the solid Cc!) 


Conduction parameter ( ee 
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BTU/hr°F ft 


BTU/1b°F 


inches 
inches 4,0 
lbs/£t> 


inches Hg 


lb/hr ft 
lb/hr ft 
sec/inch 
inches 11,0 


inches H,0 
(inches) + 
°F 


language along with 


lb/hr 


BTU/sec °F 


BTU/ °F 





SLOMAX Maximum slope of generalized cooling curve - 


G Mass flow velocity 1b/hr ace 
RNUMFR Reynold's number for friction = 
FFR Fanning Friction factor - 


Reduction of data is not complete since j remains to be found, 
With SLOMAX and CONDPAR, enter Table VIII, and find Neu The Colburn j 
factor follows from equation (A-22), Appendix A. 

The program is flexible in that an unlimited number of matrices 
may be reduced by changing the first card of the program appropriately. 


(L = the number of matrices to be reduced.) 


oy 


- itt = <a. > i=» 


nmin. 4 § Mia 


2, f+ ee —-= = 


(oo € OEE he =——« o@ ¢ 
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